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We present the formalism and demonstrate the use of the overlapping 
muffin-tin approximation (OMTA). This fits a full potential to a superposi- 
tion of spherically symmetric short-ranged potential wells plus a constant. 
For one-electron potentials of this form, the standard multiple-scattering 
methods can solve Schrodingers' equation correctly to 1st order in the po- 
tential overlap. Choosing an augmented-plane-wave method as the source 
of the full potential, we illustrate the procedure for diamond-structured Si. 
First, we compare the potential in the Si-centered OMTA with the full po- 
tential, and then compare the corresponding OMTA N-th. order muffin-tin 
orbital and full-potential LAPW band structures. We find that the two 
latter agree qualitatively for a wide range of overlaps and that the valence 
bands have an rms deviation of 20 meV/electron for 30% radial overlap. 
Smaller overlaps give worse potentials and larger overlaps give larger 2nd- 
order errors of the multiple-scattering method. To further remove the mean 
error of the bands for small overlaps is simple. 

PACS numbers: 71.15.-m, 71.15.Ap 

1. Introduction 

The linear muffin-tin orbital (LMTO) method in the atomic-spheres ap- 
proximation (ASA) [TJ [5] has been used for over three decades as an ab initio 
method for electronic-structure calculations. One of the defining characteristics 
of this method is its small, localized [3j basis set which provides highly efficient 
computation and does not require crystalline symmetry [4]. Equally important 
is the fact that the similarity of LMTOs to atomic orbitals makes the method 
physically transparent. This means that it is easy to extract e.g. information 
about the nature of chemical bonds, the symmetry of band states, a.s.o.. 

The extreme speed and simplicity of the LMTO-ASA is to a large extent 
due to the ASA which takes the one-electron potential and charge density to be 
spherically symmetric inside space-filling Wigner-Seitz spheres whose overlap is 
neglected. This approximation is very good, not only for closely-packed solids but 
also for open systems in which the interstital sites have high symmetry and can be 
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taken as centers for " empty" Wigner-Seitz spheres containing no protons but only 
electrons. The diamond structure is such a case, because both its potential and 
charge density are are approximately spherical inside equal-sized spheres packed 
on a body-centered cubic (bcc) lattice. However, for open systems with atoms of 
different sizes and/or with low symmetry, it can be cumbersome, if not impossible 
to find a good ASA. 

Since most systems of current interest are low-dimensional and multi-component, 
and since the cost of computation has followed Moors' law for decades, the LMTO- 
ASA method is now mainly used to extract information about the electrons from 
accurate calculations which use one of the plane-wave supercell methods. These 
methods employ the variational principle with large plane-wave basis sets, either 
for the full-potential Schrodinger equation with linear augmented plane waves 
(LAPWs) [TJ [5] or projector-augmented plane waves (PAWs) [5J, or for a pseu- 
dopotential with straight plane waves [7]. In this mode of operation, one adjusts 
the LMTO-ASA sphere sizes and empty-sphere positions by trial and error such 
as to obtain the best possible agreement with the band structure obtained with 
an accurate method. 

Attempts have been made to improve the accuracy while preserving the 
desirable features of the LMTO-ASA. First [5], it was realized that the multiple- 
scattering method of Korringa, Kohn and Rostoker (KKR) [9], the mother of 
all MTO methods, solves Schrodingers' equation not only for a MT potential, 
that is a potential which is spherically symmetric inside non-overlapping spheres 
and flat (=0) in between, but also for a superposition of spherically symmetric 
potential-wells to 1st order in their overlap, with no change of the formalism. 
This explained why the ASA had worked much better than the MT approxi- 
mation (MTA) and begged for the use of even larger potential spheres. The 
overlapping muffin-tin approximation (OMTA), however, turned out not to work 
with conventional LMTOs because the energy-dependence of their envelope func- 
tions, the decaying Hankcl functions, had not been linearized away in the same 
way as the energy-dependence of the augmenting partial waves inside the spheres, 
but put to a constant (e.g. =0). As a consequence, exact MTOs (EMTOs) with 
envelopes of energy-dependent, screened solutions of the wave equation had to be 
found, and then linearized [10]. With the resulting, so-called 3rd-generation LM- 
TOs (LMT03s), it became possible to solve Schrodingers' equation for diamond- 
structured silicon using the OMTA without empty spheres [TT]. Specifically, the 
ASA potential calculated self-consistently with empty spheres (and the 14% ra- 
dial overlap of bcc WS spheres) was least-squares' fitted to a superposition of 
potential wells centered only on the Si sites, whereafter Schrodingers' equation 
was solved for this OMTA potential with the LMT03 method. By varying the 
OMTA overlap, a minimum rms error of 80 meV per valence-band electron was 
found around 30% overlap: for smaller overlaps the OMTA-fit was worse, and 
for larger overlaps the 2nd-order overlap error of the LMT03 method dominated. 
This ability of the LMT03 method to compute the band structure for an OMTA 
potential was not developed into a charge-sclfconsistcnt method, despite several 
attempts |T2J |T3] , but for the EMTO Greens '-function method this was done by 



3 



Vitos et al. who devised a spherical-cell approximation for closely-packed alloys 
and demonstrated its efficiency for overlaps around 30% |14j . Instead, the unique 
ability of the LMT03 method to generate super-minimal basis sets through down- 
folding was explored, but that led to the next problem: The smaller the basis set, 
the longer the range of each MTO, and the stronger its energy dependence. For 
very small basis sets, this energy dependence may be so strong that linearization 
is insufficient. That problem was finally solved through the development of a 
polynomial approximation of general order (N) in the MTO Hilbert space, and 
the resulting JVMTO method [15] has now been in use for nearly ten years, in 
particular for generating few-orbital, low-energy Hamiltonians in studies of cor- 
related electron systems, but always in connection with potentials generated by 
the LMTO-ASA. 

Here, we shall present and apply the least-squares' formalism to fit a full po- 
tential in the form delivered by any augmented plane-wave method to the OMTA 
form. This procedure is demonstrated for diamond-structured silicon, first by 
comparison of the full with the OMTA potential, and then by comparison of the 
full-potential LAPW and OMTA 7VMTO band structures. Our current strategy 
for extracting electronic information from an accurate plane-wave calculation, is 
thus to fit its output potential to the OMTA and then perform 7VMTO calcula- 
tions for the OMTA potential. This procedure is more simple and accurate than 
going via the LMTO-ASA. 

2. The Overlapping Muffin-Tin Approximation 

In this section we shall explain the least-squares' procedure for obtaining the 
OMTA to a full potential. In the first two subsections we shall present the general 
formulas 13J, and in the last we shall present results for the specific analytic form 
of the full potential provided by an augmented plane-wave method. 

We wish to minimize the value of mean squared deviation, 
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of the the full potential, F (r) , from its OMTA counterpart, g + J^j fj( r j)- The 
spherically symmetric potential wells fj(rj) vanish for rj > Sj and are centered 
at the sites Rj of the atoms, fl is the (cell) volume. To be determined are thus 
the radial functions fj (r) and the constant g, for given F(r), £1, Rj, and Sj. 

2.1. Determining the spherical potential wells, fj(r) 

Variation of AF with respect to fj (r) leads to the equation: 







0( Sj - r) (/• •••'• -g)- fj(r) - 6( Sj - r) £ . jf (2) 



J 

centered at site j and of radius r. Eq. (|2|) thus translates into the requirement that 



where F J,r and fjf are the averages of the respective functions over the sphere 
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the spherical averages of F(r) — g and Ylf fj' ( r j' ) around site j be the same. The 
problem of calculating F J ' r for the LAPW potential will be addressed in Sec. 12.31 
By reordering of the terms in Eq. ([2]), and evaluating the spherical averages, we 
obtain: 



rfj(r) = 9{ Sj - r) 



(F- -g)-Y, ^ + ^~ d ^ f" r'f f (r')dr' 

za jj' Jd H t-r 



j'^j 

"(3) 

where djji = |Rj — Rf\ is the distance between respective sites. We have assumed 
that Sji < djj' , i.e. that a potential sphere never includes the center of another 
sphere (moderate overlap). Repeating the procedure for all inequivalent sites, 
Eq.s ([3]) become a set of linear integral equations which may be solved for a given 
g = go by self-consistent iteration. The iteration can be initiated e.g. by taking 
fj,(r') = 6( Sj , - r') (P'' r ' - g) on the right-hand side of Eq. ©. 



2.2. Determining the constant background, g 

A similar expression is obtained for the optimal value of the constant g, 
namely: 

0-§- g AF = F-g-J2Jl (4) 
where F and fj are the volume averages of the respective functions: 

P = h L F(r)dr ' and Ij = ^I a /i(|r ~ Rjl)dr = ^ £ f ^ r)r2dr (5) 

Just like before, Eq. (j4j expresses the condition of equality of averages of the full 
and OMTA potentials. 

We need to solve the coupled f and g equations. Since these are linear, we 
may use the superposition principle to obtain: 

fj (5; r) = fj (g ; r)-(g- g ) hj (r) , 

where the " structure functions" hj (r) are the solutions of the /-equations ([3]) for 
F (r) — g = 1. The optimal g-value is thus: 

9 = 90 + (P - go + U (9o)) /(l - J2j h 

2.3. Averages of the LAPW full potential 

The procedure outlined in section 12.11 requires calculation of the volume 
average F, and the spherical averages, F J ' r , for r < Sj < min^/ djji around all 
inequivalent sites of the full potential. This task can be achieved with relative 
ease thanks to the analytical form of the potential delivered by an augmented 
plane-wave method: Space is divided into a set of atom-centered, non-overlapping 
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augmentation spheres with radii, a,-, and the interstitial. Inside an augmentation 
sphere the full potential is given by the spherical- harmonics (Y/ m ) expansion: 



F(r) = ^2F jlm (r j )Y lm (r j ) for rj < 



(6) 



where i\ 



R, 7 " 



r, rj. In the interstitial, the full potential is described as a 



linear combination of plane waves which we shall let extend throughout space: 
F(r) = Yl 6 ( a i - r i) E ( F 3i™(rj) - E r F * lm (G) Ylm ^ 



lm 
Gr 



(7) 



The first term here provides the augmentation of the plane waves, i.e. it ads 
expression ^ and subtracts the corresponding plane- wave part inside the spheres. 
G = GG are the reciprocal lattice vectors and F (— G) = F* (G) because F (r) 
is real. Moreover, j; are the spherical Bessel functions, and 



F jlm (G) ee W ]>\ F (G) e iG ^y^(G). 



(8) 



The spherical average around site j is the spherical component of the po- 
tential ([6]) inside the augmentation sphere: F J,r = Fjqq {r) /V™, if r < a,j. If 
7' > a-,- and the r-sphere at site j cuts into the augmentation-sphere at site j', 
the spherical average of the corresponding term in the first line of Eq. ([7]) is most 
easily evaluated when, for the spherical-harmonics expansion inside the j'-sphere, 
we turn the z'-axis to point towards Rj, i.e. along the djj, direction: 

Y lm (r r ) = J2 m , D l mm , (djj,) Y lm , (fj, • djj,) ee J2 m/ D l mm , (djj,) Y lm ,{# ,<p% 

because then the integral over <p' selects the term with Y/o (6', ip') = \J~^^Pi (cos 9') 
and only the integral over 9' remains. Including now also the spherical average of 
the second line in Eq. 0, the result becomes: 



p3- 



r'F jfl (r')Pi 



r' 2 — r 2 + djj, 



+ ^io (^r)Fjoo (G)/V4^, 



efr' 



(9) 



where jo 0*0 — sinx/x and 

121 + 1 



F ]3 , l {r')^9{a 3 



= 9 (aj, - r') 
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co=30% 



w=60% 




along [100] 



along [100] 



along [100] 



Fig. 1. Full potential for diamond-structured Si along the main crystallographic axes 
(dashed) and its OMTAs for 10, 30, and 60% radial overlaps (full). The double-headed 
arrow indicates the range of the valence bands. 



Here, Fjj, w (G) is given by Eq. © with Y* {G) substituted by 
Finally, the volume average of the full potential 10 is 



21+1 P 
in rl 



Gd 



J 3 




Fjoo{r)r 2 dr 



"3 A^i 

G 



+ F(G=0) ; 



where j\ (x) — (sinx — xcosx) /x 2 . 

3. Example: Diamond-structured silicon 

In this section we shall compare the potential and band structure calcu- 
lated self-consistently using the full potential LAPW method as implemented in 
WIEN2k iTB] with respectively the OMTA to this full potential and its band struc- 
ture calculated using the iVMTO-OMTA method. As in the previous (AE)SA-i;s- 
OMTA LMT03 study [11 , in the present full-ws-OMTA NMTO study we choose 
diamond-structured Si and OMTAs with atom-centered spherical wells in order 
that the study be relevant for open structures in general. 

Fig. [T] compares the full potential with its OMTA for radial overlaps, u> = 
(2s — d)/d, of 10, 30, and 60%. We see a marked improvement in the quality of 
the fit as the overlap, i.e. the range of / (r) , increases. With 60% overlap the 
OMTA follows the full potential not only inside the potential wells, i.e. close to 
the nuclei where the spherical approximation is at its best, but also along the 
bond between the Si atoms and in the interstitial space. 
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Fig. 2. Rms errors of the iVMTO-OMTA valence and four lowest conduction bands as 
functions of the radial overlap, with respect to the full-potential LAPW bands. 

From this, one might believe that it is always preferential to use large 
overlaps. But this is not so, because multiple-scattering methods are only "self- 
correcting" to lst-order in the potential overlap (oc uj 2 ) [H1[TT]. For large overlaps, 
the errors of higher order (oc oj 4 and higher) outweigh the effect of the improved 
description of the potential. We thus expect that there exists an optimal value 
of the overlap which yields the best accuracy of the calculated band states. In 
Fig. 2 this is confirmed by explicit calculations of root mean squared (rms) error 
of the 7VMTO-OMTA bands with respect to their full-potential LAPW counter- 
parts as a function of the radial overlap. The two curves are for respectively the 
four valence and the four lowest conduction bands. For both types of bands, the 
error initially decreases as a result of the improved OMTA. For the valence bands 
the optimum of 3 mRy/band = 20 meV/electron is reached at 30% overlap, after 
which the error increases sharply, whereas for the conduction bands the optimum 
of 50 meV/e is reached at 40%, wherafter the multiple-scattering errors increase 
slowly. This difference in behaviour is due to the valence states being bonding and 
the conduction states antibonding, whereby the former probe the overlap region 
more than the latter. Whereas the rms-error curve for the conduction band is 
nicely flat, that of the valence band is not. 

Fig.[3]shows the full-potential LAPW band structure and the A^MTO band 
structures for the three OMTA potentials in Fig. [TJ Excellent agreement is 
achieved with 30% overlap for both conduction and valence bands. For smaller 
overlaps, the iVMTO-OMTA bands lie to high, the valence bands in particular, 
and for larger overlaps they lie too low. The first point follows from the varia- 
tional principle and the second from the fact that the multiple-scattering error 
oc uj 4 is negative definite. This was explained in the previous study [TT] , which also 
demonstrated the use of various corrections. The simplest of these is to remove the 
mean error of the energies for a group of bands, say the valence bands, by weight- 
ing the squared deviation in Eq.Q] with the charge density of those bands, and 
then keep the weighting in the g-equation only (see Sect. 12. 2|) . Since the LAPW 
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Fig. 3. Band structure of Si calculated with the full potential LAPW method (dashed) 
and with the AMTO-OMTA method for 10, 30, and 60% radial overlap (full). The 
iVMTO basis set consisted of the 9 s, p, and d Si-iVMTOs with the fs downfolded and 
had N—2 with expansion energies -0.5, 0.3, and 0.6 Ry [15] . 

valence charge density has the same form as the full potential, implementation 
of this requires use of formulas given in Sect. 12.31 This correction will move the 
iVMTO-OMTA bands down by the appropriate amount for small overlaps and, 
hence, considerably lower the rms-error curve for small overlaps. An elegant cor- 
rection of the large-overlap errors remains to be found; the straight-forward but 
clumsy one is to evaluate the kinetic energy properly in the overlap region for the 
ATViTO basis. 
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